Gauge invariance and electron spectral functions in underdoped cuprates 
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The single particle spectral function for the normal state of underdoped high T c cuprates is 
studied within the slave particle framework. We find that the presence of a massless dynamical 
gauge field - a direct consequence of the quantum order - explains the broad, but not totally 
incoherent, line-shapes observed in experiments. The issue of the negative anomalous dimension of 
a recently proposed gauge invariant single particle amplitude is also considered. We show how the 
anomalous behavior of the single particle amplitude can be incorporated within the slave particle 
approach and, thus reinterpreted, lead to physical phenomenology. 
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I. INTRODUCTION 

In recent years gauge theories have become a promi- 
nent tool in the description of strongly correlated electron 
systems |lj g, || . The gauge field appears because we at- 
tempt to use fermion operators to describe excitations of 
a strongly correlated state. 

We know that the fermion operators are useful at low 
energies only when the ground state wave function can 
be described by a Slater determinant. If the ground 
state is a strongly correlated state, we cannot use the 
fcrmion/boson operators to write down a free low en- 
ergy effective theory. However if we insist on writing ev- 
erything in terms of fermion/boson operators, one finds 
that we can use fermion/boson operators to describe the 
low energy physics of a strongly correlated state provided 
that gauge fields are introduced. This is why low energy 
effective theories of strongly correlated states, such as the 
quantum Hall states and spin liquid states, are all gauge 
theories. 

Realizing that the electron operator is not a good start- 
ing point to describe a strongly correlated state, we can 
rewrite the electron operator as a product of several other 
operators in order to study those states. These operators 
are called parton operators (the spinon and holon oper- 
ators are examples of parton operators). We then con- 
struct mean- field states in the Hilbert space of partons. 
After identifying the gauge structure - the transforma- 
tions between partons that leave the electron operator 
unchanged, we can project the mean-field state onto the 
physical Hilbert space (i.e. the gauge invariant space) 
and obtain a strongly correlated electron state. Although 
the parton mean-field state is a Slater determinant, the 
projected state is not and can be used to describe the 
correlations in certain strongly correlated states. This 
procedure, in its general form, is called projective con- 
struction, which is a generalization of the slave-boson 
approach^ || [| |j. The generic projective construction 
and the associated gauge structure was discussed in de- 
tail for quantum Hall states in j?J . The resulting effective 
theory arising from the projective construction naturally 
contains a gauge structure. In this paper, we will use the 
projective-construction/ slave-boson-approach to study 



the electron spectral function in underdoped high T c su- 
perconductors. 

One convenient way to describe dynamical gauge theo- 
ries is via the path integral approach. Traditionally, how- 
ever, the formalism which allows for most direct compar- 
ison with experimental results is based on Green's func- 
tions and diagrammatic techniques. This turns out to be 
a kind of nuisance since the standard expression for the 
single particle Green's function (the cornerstone of dia- 
grammatic perturbation theory) is not gauge invariant. 
The lack of gauge invariance makes it hard to extract any 
meaningful quantitative results based on these Green's 
functions. 

In the first part of the paper we present a calcula- 
tion of a particular gauge invariant Green's function for 
2+1 dimensional massless Dirac particles coupled to a 
dynamical gauge field. The motivation for this analy- 
sis arises from an effective continuum description - the 
Algebraic Spin Liquid (ASL)|||, |{|, which we believe un- 
derlies the strange phenomenology of the normal state of 
underdoped high temperature superconductors. Within 
this continuum model, which is rooted in the slave-boson 
approach to strongly correlated systems, the spin de- 
grees of freedom are carried by massless Dirac particles 
called spinons. Our early attempt to interpret the re- 
sulting spinon propagator as describing the propagation 
of the physical particle is however thwarted by the ap- 
pearance of a negative anomalous dimension for the par- 
ticular gauge invariant spinon amplitude at the ASL fixed 
point. On reinterpreting the proposed amplitude as rep- 
resenting a two-particle spinon-holon propagator, we find 
that the main effect of the gauge fluctuation is to bind 
the two particles (via its effect on the vertex) to a more 
coherent entity. Thus we are led to the heuristic descrip- 
tion of the physical electron propagator presented in the 
last section of the paper where in addition to the spinon 
contribution we include the bosonic holon contribution. 
Not surprisingly the resulting spectrum has no quasipar- 
ticle peak, however it is also not totally incoherent once 
the gauge fluctuations are included. The emerging pic- 
ture for the electron propagation is thus one of spinons 
and holons (separate at the mean field level) whose at- 
tractive interaction mediated via the gauge field leads to 
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the formation of a more coherent structure on top of the 
incoherent mean-field background [fl()| . In analyzing this 
problem we have been guided by a first quantized path 
integral approach to the physical Green's function which 
makes the gauge dependence of the amplitude particu- 
larly transparent. 



II. PATH INTEGRAL FORMULATION 

Before plunging into the second quantized formalism 
let us look at the problem in the light of single particle 
quantum mechanics. Our starting point is the following 
definition for a gauge-invariant Green's function in first 
quantized notation 

G(x) = ( DaDxe-^ dT{L[x(T) - a]+L[a]) e- lS \ a » dx » (1) 



where | is the straight line connecting to x, fi = 0..2 
and the metric is Euclidean. The action for the gauge 
field is given by 



S\a,, 



which is obtained by integrating out massless Dirac 
fermions (the spinons) [See Eq. (ph]. The action for the 
particle moving in an external gauge field can be taken 
to have the following form 



d,TL[x(T) 7 a}= / <1tL[x{t) 
Jo 



ci^dx ^ 



(3) 



where J indicates a path from to x described by x{t). 
L[x{t)\ remains unspecified since we will be concerned 
with different types of particles. Putting the above to- 
gether we obtain the following expression for the Green's 
function 

G(x ) = / DaDxe~ dT{L[x ^ +L ^ ) e i ^ a * dx *~ i ^ a * dx > Jl {^) 



This expression exhibits the gauge dependence in a par- 
ticularly appealing form where the line integrals close 
up into a contour, of which one side is the propagat- 
ing path of the particle and other side the straight-line 
path needed to make the gauge invariancc manifest. Via 
Stokes theorem we see that the sum of the line integrals 
is nothing but the flux going through the contour. So far 
our expressions have full rotational symmetry in order 
to perform the calculation we now chose to parameter- 
ize the particle world-lines by their xq = r coordinate. 
Hence J a^dx^ goes over to J (a(x)^ + ao)dr where a, x 
denote the 2 spatial components of the corresponding 
3-vectors. In order to simplify our analysis further we 
consider the propagation of the particle only along the 
time direction i.e. G(0,r). Considering the propagation 
between equal spatial coordinates we have no contribu- 
tion from the 2-velocity component along the straight 



line path which simplifies this part of the phase factor 
to J, aodr. Next we want to perform the average over 
the gauge fluctuations which is straightforward given the 
Gaussian weight (0). Reading off the gauge propagator 



D (a^a v (k) 



(2tt)-- 



q 2 



In order to invert the polarization we chose the ao = 
gauge which yields 



D 



(a^q)a u {k)) = {2nf5{q+k)5^ [S ZJ + ^ 



where i, j — 1,2 denote the spatial components of q. Note 
particularly that in choosing the ao = gauge we got rid 
of the remaining contribution of the straight line path to 
the phase factor in the expression for the Green's func- 
tion. Now we perform the average over the gauge fluctu- 
ations 



(exp(- 



axdr)) £ 



exp 



x(r)x' '(T')dTdT' e m e ik *' (a(g)a(jfc)), 

q.k 



The calculation is straightforward albeit some care has 
to be taken with the boundary conditions. The averaging 
results in 



x(r)x'(r') + l 



x- 



+ (t-t') : 



(t — T')dTdr' x 
1 

x 2 + (r - r') 2 
1 



+ (r-r') 2 {t-t'Y 



(5) 



The last three terms in the above expression (^|) arise 
from the boundary conditions and thus explicitly depend 
on the direction of propagation. In order to extract the 
dependence on Xf — xi explicitly we can now exploit the 
rotational invariance in the Euclidean formulation to re- 
express the above in vector notation. 




exp 



dx ■ ndx' ■ n 



Xi-ri J Xi-n 

dx ■ dx' 



'(x • n — x ■ ft) 
dx ■ ndx' ■ ft 



(x — x') 2 [x — (x' ■ ft)n\ 2 
dx ■ ndx' ■ ft 



[x' — (x ■ ft)n\ 2 [{x ■ n)n — (x' ■ ft)n] 2 



(6) 



with 



= the unit vector along the classical 

straight line path. Expressions (||) and (^|) exhibit nicely 
the retarded nature of the effective interaction. In prin- 
ciple, to obtain the propagator G(x), we would now have 
to perform the average over the trajectories with the free 



particle action appropriate for the particle in question. 
For a massive particle we can neglect the spatial part 
of the retardation since it is down by a factor of the 
speed of "light" (the spinon velocity) which we have set 
equal to one. Hence we can see that only the term pro- 
portional to x(f)x'(f') gives a non- vanishing contribution 
whose effect is to renormalize the mass of the particle. 
For a massless relativistic particle however both spatial 
and temporal retardation are equally important and the 
argument above cannot be applied. We find that the in- 
tegrals in (|B|) are all dimensionless. In the large r limit, 
those integrals have the form cif + In f- , where In is 
the short distance cut-off scale. The c\ term corresponds 
to mass generation. Within second quantization, as we 
will discuss next, the mass term cannot be generated if 
we regularize our theory at short distances in a way con- 
sistent with the underlying lattice symmetries + gauge 
structure ]TT[ ] . Thus we can set the regularization depen- 
dent term C\ to zero. The true leading term has a form 
c^lnr (02 is regularization independent). Those contri- 
butions can only change the exponent of the free Green's 
function: 1/t 2 — > 1/t( 2 ~ q ). The effect of the gauge inter- 
action is to modify the exponent in the algebraic decay 
of the Green's function, which means that the gauge in- 
teraction is a marginal perturbation. 



III. MASSLESS DIRAC SPINORS 

Having used the first quantized version of the path in- 
tegral to guide our intuition we will henceforth return 
to field theory and in particular to the problem of mass- 
less Dirac spinors coupled to a U(l) gauge field. The 
Euclidean model which we are concerned with is the fol- 
lowing 



= / DfD^expi - 



N 
cr=l 



(7) 



with the dynamics for the gauge field given by 

Z a = J Da^cxp^-i J -^^a^{q)\\^a u {-q) 



n, 



-w^ <w 



q 2 



(8) 



The Fermi field $ ff is a 4 x 1 spinor = (^nVva) 
with ipv i = ({ <Tle ) ,Vv.2 = (l° 2 ")- The notation (l,2,e,o) 
arises naturally when the theory (|7|) is derived as the low 
energy effective description of the so called 7r-flux phase 
of Hubbard type lattice models fl2fl . There 1,2 denotes 
two types of Fermi points and e,o denote even and odd 
lattice sites respectively. Notice that we have generalized 
to N fermion species in order to have a more controlled 
perturbation theory (N=2 in the physical case). The 4x4 
7^ matrices form a representation of the Dirac algebra 



E(p) .f s q 
_ ^ * >. _ 

— 

- - — <— 



p-q 

FIG. 1: Exchange diagram used in the calculation of the self- 
energy 



{7/i, 7i/} = 2(5 M i/ (/i, v = 0, 1, 2) and are taken to be 

7.= (S ° Y 7, = f? _° A (9) 



72 



o-i 

-01 



with ct m the Pauli matrices. Finally note that = 

^From the path integral discussion above we have 
learned a way to calculate a particular gauge invariant 
Green's function for the Fermi fields. Note that a certain 
choice for the gauge has to be made before we can invert 
the polarization operator to obtain the gauge propagator. 
In order to get rid of the contribution from the straight 
line path in G, we can fix the gauge in such a way as 
to make it vanish along that very line. In this particu- 
lar gauge, G has a form of a normal one-particle Green's 
function and can be calculated using the usual method 
through the self energy: 

G(p) = G (p) + G (p) (-£($) G (p) + ■■■ 

where £ to lowest order in 1 /N is the exchange diagram 
depicted in Fig.Q. However, the above trick only works 
for one particular direction x in real space. For differ- 
ent directions, we need to choose different gauges. We 
can put this into a more precise mathematical form. By 
choosing a gauge h ■ a = 0, we can calculate the self en- 
ergy and the corresponding Green's function (which is 
denoted as Gn(xj) through diagrams. However, Gn(x) is 
not equal to the gauge invariant Green's function G(x). 
But the two Green's functions are closely related: 



G(x) — Gfi(x)\n = x/\S\ 



(10) 



In the following, we will calculate G through calculating 
Gfi—f • 

To leading order in the gauge field fluctuations the 
fermion self-energy is 

where Go = ~'^2 7 is the free fermion Green's function. 
The sign convention for the self energy has been chosen 
such that 

G- l {k) =tVy* + £(A) 



4 



To invert the gauge polarization, we will choose the 
a = gauge. This implies that our calculation gives a 
correct result for the equal-space propagation only. Em- 
ploying relativistic invariance however, we can immedi- 
ately extend the result to all of real space. With the 
above gauge choice 



8 1 



s^Suj [Sij + -V- 



and we arrive at the following expression for the self- 
energy at T=0 



(ii) 



+ Y 



+ Y 



d 3 q 
(2tt) 3 

d 3 q 
(2^3 

d 3 q 
(2tt)3 



(?o -Po)q 2 + 2(g -Paha 



(Px - qx)(q 2 x - ql) + 2q x q y (p y - q y ) 



V^ip- q) 2 ql 

(Py ~ %){<& ~ ll) + 2qxq y (Px - q x ) 



<f (p- q) 2 q 2 



In order to evaluate the self energy we need to take care 
of both the UV divergences usually encountered in loop 
integrations and the singularities arising from the l/<?o 
terms which are a hallmark of the incomplete gauge fixing 
in the temporal axial gauge. Details of this calculation 
are presented in the appendix. We find for the self-energy 



v 16 1 o; ( A 



(12) 
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Had the self-energy taken the form 



A 



-i^p j°ln 



(13) 



A 

m 

X; ; i/7) - 'In ( A 



we could read off the anomalous dimension of the fermion 
field and found the Green's function Gf (x) 



Gf(x) ex 



(14) 



£3+27* 

Since the above Gf has no explicit f dependence, we find 



G(x) = G f {x) cx 



£3+27* 



(15) 



However, the coefficients of the two terms in Eq. ( |l2] ) 
are not equal and the above calculation of G(x) is not 
valid. But noting that the coefficients in Eq. ( |l2| ) have 




P-0 



p_0 - q_0 



FIG. 2: Exchange diagram used in the calculation of the bo- 
son self-energy 



the same sign - corresponding to a negative 7^ 
allows us to conclude that Gf (r) has the form 



this 



Gf(r) cx 



T7 



16 1 



1 



-3+27* ' 



-iV3^- 7 *--iV3^ (16) 



We see that the effect of gauge interactions is to make the 
propagator decay slower. In the path integral picture, the 
gauge interactions bring the particle path closer to the 
straight-line path so the area enclosed by the two paths is 
smaller. Since the particle path is closer to the straight- 
line path in the r direction, this leads us to believe that 
the second pi term in Eq. (n3) is not important compared 
to the first po term. Hence 



16 1 

= — «TT 

N 6TT Z 



(17) 



This gives us 



G 



2-3+27* ' 



16 1 



^ = -N3^ (18) 



We would like to mention that the renormalization pro- 
gram within the temporal axial gauge is not established 
(see Ref. |f| for a discussion on renormalization within 
non-covariant gauges) . To check and to confirm the result 
Eq. ([l8]), we follow a different route where we consider 
in addition to the massless fermions coupled to the gauge 
field a heavy boson field whose "current" will take care of 
the straight line path in the equal space correlater. This 
type of formulation is borrowed from heavy quark effec- 
tive theory approaches to QCD Jl4j]. The task at hand 
is thus a calculation of the wavefunction renormalization 
of the boson fermion bilinear j = tyb with the following 
Lagrangian 



{dyrf + ia^)^ a + b* (do + ia Q )b + C a 



In order to calculate the wavefunction renormalization 
of the above composite we need the field strength renor- 
malizations Z^ , Zb for both the fermion and the boson 
field as well as the vertex renormalizations Zr- To obtain 
the boson wavefunction renormalization we calculate its 
self energy due to the gauge field fluctuations FigQ. As 
can be seen from the above Lagrangian, the boson being 
static only couples to the Coulomb field. Thus we arrive 
at the following expression for the boson self energy. 



^b(pa 



d 3 q ijoD o(q)ijo 
(2tt) 2 i(p Q - g ) 



(19) 
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FIG. 3: Boson- Fermion vertex with the po flowing along the 
static boson line 



whose divergent part evaluates (see appendix) to 

1 



zE 6 (p ) = Pa 



ir 2 N 3-d 



(20) 



where the sign convention for the self energy has been 
chosen in accord with the one for the fermion field men- 
tioned above. 

Thus we find for the wavefunction renormalization the 
minimal 



Gb(pa) 



ip a + S(p ) 



Zb_ 

ipo 



1 



1 



TT 2 NZ-d 



This result has been obtained in the Landau gauge. The 
fermion wavefunction renormalization can be obtained 
similarly from the self energy and reads in the Landau 
gauge 



= 1 



1 



3ir 2 N 3 - d 



(21) 



Finally we need to obtain the vertex renormalization. 
From Fig. II we can read off the following expression for 
the vertex correction 



d 3 l i<y»(-i)'f(k + l) e 18 0u iD lu ,(!) 



(2tt) 3 



(k + l) 2 i{po + lo) 



where 1 is the fermion-boson vertex whose renormaliza- 
tion is sought. By simple powercounting we see that 
the loop integration is only logarithmically divergent and 
hence we can evaluate it by setting all external moment 
equal to zero which leaves us with 



1 



d 3 l 
(>]3 

= 1 



7o7% _ l^j%J>xV 
l 3 lo 
dH 



(2tt) 3 



1 



Pi 

I 5 



where the first term was simplified by the fact that only 
the 7°/o term survives integration over spatial momenta. 
Thus we conclude that within the Landau gauge the ver- 
tex doesn't renormalize to one loop order. 



We now have all the ingredients to calculate the wave 
function renormalization of the fermion boson bilinear in 
the minimal form 



Zj = zl ,2 z\l 2 z T = 1 



16 1 



1 



N 3tt 2 3 - d 



(22) 



This result is gauge invariant. From the above discussion 
we can see that within the axial temporal gauge there 
is no Coulombic interaction and hence both the vertex 
correction and the boson wavefunction renormalization 
are absent which would suggest that 



Zv(a = 0) = Z 2 = l 



16 1 



N 3tt 2 3 - d 



which yields 



7* = 77 



16 1 

N 3^2 



(23) 



for the anomalous dimension of the fermion field. 

An important point to note here is the fact that the 
gauge fluctuations do not generate a mass or a chemical- 
potential term perturbatively. In a diagram calculation, 
a regularization dependent mass/chemical-potential term 
actually appears in the self-energy. We have set such a 
regularization dependent mass/chemical-potential term 
equal to zero. In the next section we will see that our 
Dirac-[/(l) gauge system is the low energy effective de- 
scription of the so called staggered flux (sF) phase (the 
mean-field phase corresponding to the pseudogap regime 
within the SU(2) slave boson theory). The lattice sF 
phase provides a short distance regularization. Now the 
question is whether such a regularization generates a reg- 
ularization dependent mass/chemical-potential term or 
not. One way to answer the above question is through 
a symmetry argument. In Ref. it is shown that 

any change of the mean-field sF ansatz cannot generate 
a mass term or chemical potential term if the change does 
not break any symmetries and does not break the U(l) 
gauge structure. This implies that if the regularization 
is consistent with the underlying lattice symmetries and 
does not break the U(l) gauge structure, then the reg- 
ularization cannot generate any mass/chemical-potential 
terms. This is the reason why we can set the regular- 
ization dependent mass/chemical-potential term equal to 
zero. We see that the massless Dirac fermion is pro- 
tected by the lattice symmetries and by the U(l) gauge 
structure (or more precisely, protected by the quantum 
order[jll|). Thus what we really have learned from the di- 
agram calculation is that the mass or chemical-potential 
is not generated spontaneously from the infrared diver- 
gences (at the perturbative level). Combining the sym- 
metry consideration for the high energy cutoff and the di- 
agram calculation for the infrared, we find that no mass 
or chemical-potential is generated from the self energy 
corrections to the fermion Green's function. In contrast, 
the first quantized path integral analysis does not know 
about the underlying lattice symmetries and, as a result 
allows for corrections to the mass or chemical-potential 



6 



due to dressing with the short-distance gauge fluctuations 
(see discussion in section ||). We would like to mention 
that the Dirac-£/(l) gauge system was studied in detail 
in Ref. M| using a different gauge fixing. It was shown 
that the mass/chemical-potential term is not generated 
perturbatively to all orders in a 1/TV expansion. Further- 
more note that since the conserved current (that couples 
to a^) cannot have any anomalous dimension, the gauge 
interaction is an exact marginal perturbation. 

Hence we can take the anomalous dimension seriously 
and obtain the full Green's function in the form 



G(fc) 



-iC 



£2-27,1, 



TV 37r z 



where C is determined by the energy range over which 
the effective theory is supposed to be valid. Comparing 
this dressed propagator with the free fermion Green's 
function Go = ~' fc 2 7 we see that the inclusion of the 
gauge fluctuations has destroyed the coherent quasipar- 
ticle pole. 

At this point a couple of comments are in order. So 
far we have not addressed the fact that the sign of the 
anomalous dimension (leading to a slower decay of the 
single particle propagator in real space) is highly counter 
intuitive. Fermion-fcrmion interaction is normally ex- 
pected to result in a suppression of the single particle 
propagator - not an enhancement. This fact combined 
with wishful thinking has led to the QED3 picture of the 
single particle spectrum proposed in Ref. || . This paper 
initiated a debate about the true value and sign of the 
gauge invariant anomalous dimension in massless QED3 
which is ongoing |l6|, [l?]] . The resolution of the question 
of the right gauge invariant "dressed charge" propaga- 
tor (as Khveshchenko puts it |Ti| ) is an interesting and 
important question for QED3 and beyond the scope of 
the current analysis. However, what we must conclude 
from the above stated calculation is the fact that the 
gauge invariant object which was studied here - a mass- 
less Dirac fermion coupled to an infinitely heavy particle 
(the straight-line path above) at the ASL fixed point has 
a negative anomalous dimension. In hindsight this can 
be understood as a sign of confinement - the effect of the 
gauge fluctuations is the tendency to form a bound state 
between the two particles. This picture of confinement is 
one of central importance to the question of spin charge 
separation. 

In the following we will outline a calculation of the 
electron spectral function which is more directly rooted 
in the SU(2) approach of Wen and Lee 0]. Here we in- 
clude the effect of the gauge fluctuations - the confining 
tendency, gleaned from the above analysis - more heuris- 
tically. 



IV. CONFINEMENT BY HAND 

It was mentioned above that the problem of massless 
Dirac particles coupled to a gauge field is intimately re- 



lated to a possible model description of the copper oxide 
planes in the high T c superconductors To keep the 
discussion reasonably self-contained we will first relate 
some of the details of the mean field flux state as it ap- 
pears in the SU(2) slave boson formulation M. 

The slave boson theories approach the copper oxide 
problem from the insulating side; restricting the single 
particle Hilbert space on each Copper site to three states 
with the doubly occupied one prohibited via the use of a 
Lagrangian multiplier field. Within a path integral im- 
plementation of the above mentioned constraint it is the 
Lagrange multiplier that gets promoted to the temporal 
component of a gauge field. The starting point for the 
investigation is the t-J Hamiltonian which is believed to 
contain the essential physics and exhibits nicely the com- 
petition of kinetic versus local spin fluctuations. Within 
the SU(2) formulation the physical electron operator is 
represented as follows: 



(25) 



4^ = 4(61^-4/).) 



V2 



V2 



where the following SU(2) doublets were introduced 



The ipiji,ipii are the two fermion fields representing the 
destruction of a spin up and spin down on site i respec- 
tively and hi is the doublet of bosonic fields keeping track 
of the doped holes. The procedure is now to put this rep- 
resentation into the t-J Hamiltonian 

H = PY, [ j & ■ Sj - \ n i n i) - ^U^i + h - c -)] p 

m 

which on performing a Hubbard-Stratonovich transfor- 
mation to the appropriate bosonic bond variables yields 
the following partition function 

Z = J DhDh 1 DipDd a DUe" f° >L (26) 

L = \ E Tr Ph u v] + \ E + JUv)^ 

<ij> i,i,CT 

+ E a U^yvv* + /4^) 



^2 h \ {(9t - fJ-)5ij + iUij)hj 



(27) 



The Sq fluctuations incorporate the projection to the 
space of SU(2) singlets. Furthermore note that J = 
3 J/8, t — t/2 [fil| and the matrix Uij in the form 
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contains the Hubbard-Stratonovich fields which classify 
the part of the phase diagram we are looking at. The 
mean-field phase diagram (see Introduction) is found by 
minimizing the free energy for a given number of particles 
with respect to the bond variables Uij . The phase which 
will be of interest to us is the so called staggered flux 
(sF)phase which can be represented as 



(28) 



The slave boson representation results in the following 
expression for the physical electron Green's function 

G(r» = — (T T (cj (v, t)c| (0, 0)) ) 

= -l{T r ((h^f,T)Mr,r)^\(0,0)h(0,0))) 

In the sF state we don't have any anomalous fermion- 
fermion pairing (both r 3 and 1 are diagonal in isospin 
space) and hence this simplifies to 



G(r, t) = - ~ ( (b\ (r, T )h (0, 0)/ T (r, r)/J (0, 0)) 



1 

2 

(4(r,r)6 2 (0,0)/f(r,r)/ i (0,0)) 



The angle brackets (. . .) mean averaging with respect to 
the partition function (|26|). The complication with this 
task however is the presence of light bosons whose ten- 
dency to condense prohibits sensible perturbation theory. 
Our strategy here will be to first calculate the following 
mean field decomposition 



1 



MF 



2 (6i(0,0)6 t 1 (r' ) T))(/ T (r 1 T)/T(0 ) 0)) 
I(hnfri nu t <v -^\/ft/ 



ft/ 



H (62(0,0)6T(r ) r))(/|(r' ) T)/ i (0 ) 0)) 



and then incorporate the confining effect of the gauge 
field by making the algebraic decay of the spinon Green's 
function in real space slower. The motivation for this is 
rooted in the result of the previous section (the nega- 
tive anomalous dimension for the spinon amplitude at 
the ASL fixed point) and we will discuss the implied ap- 
proximations down below. With the help of 



G n $,v n ) = £</itfy„)/l(-£,-w»)> 

G bl (k,u m ) = ^{b^q' ,v' m )b\(-k,-uj m )) 

l'> u 'rn 

G b2 (k,uj m ) = ^{b 2 (q' ,v[ n )b\(-k,-u: m )) 

9 ' ,v' m 



we can express the Fourier transform of G(f, t)mf m the 
form 



G(k,uJ n ) 



dqdv n 



(2tt) 3 
dqdv n 



Gfi (q, v n )G bl (k - q, uj n - v n ) 
3 G fi (q, v n )G b2 (k - q, u) n - v n ) 



To make contact with our previous discussions let us con- 
centrate on the fermionic part of the theory for the mo- 
ment. From the form (28) of the mean field gauge we see 
that translational invariance is explicitly broken. After 
Fourier transformation the fermionic Lagrangian is given 
in the / basis by 



E '{fMJ$(q + Q) 

\uJ n ,(T 

e/(gOcr 3 +r]f(q)a 2 



-iu> n t 



U(q) 
fa(q + Q) 



where denotes a summation over the magnetic Bril- 
louin zone which is half the size of the original one and 
takes account of the fact that in the sF phase the real 
space unit cell has been doubled with principal axis along 
±x ± y. The a matrices operate in q, q + Q space with 
Q = (tt,tt). Furthermore note that 



e f(q) = -2Jx( c os(q x a) + cos(q y a)) 
Vf(.Q) = -2JA(cos(q x a) - cos(q y a)) 
We can now simply obtain 

iv n + ef — irjf 



G n (k,v n ) 
Gfi(k,v n ) 



irjf - iv n + e/ 



Ef = e 



f 



(29) 
(30) 

Vj (31) 
(32) 



by inverting the matrix in the above Lagrangian Xa//- 
Analogously we find for the boson Green's functions 



G bl {k,uj m ) 
G b2 {k,uj m ) 



/i - e b - n] b 



El 



[iuj m - i±) 2 
iuj m - /i + e b - ir] b 



El 



4 + 



£ b{q) = -2tx(cos(q x a) + cos(q y a)) 
r)b(q) = -2tA(cos(q x a) - cos(q v a)) 

where fi is the chemical potential. 

The standard way to proceed is now to calculate the 
convolution integrals in momentum space to obtain the 
mean field Green's function. As outlined above we will 
however follow the different route of first calculating the 
Fourier transform of the above Green's functions. 

Let us first concentrate on the boson part which can 
be written as 



G bl (k,uj m ) 



' b + e b + irjb 



E b — e b - ir\ b 



2E b (iuj m - fj, + E b ) 2E b (iu m — fj, — E b ) 



(33) 
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Expanding for small k we obtain 

E b (k) ~ At X ~h{ka) 2 

e(k) ~ -4fx + ?x(H 2 

rfe(fc) ~ tAa 2 [/£-*£] 

Since — /j — i?&(0) > we can neglect the contribution 
from the first term (|33| ) and Fourier transform the second 
to obtain 



where = x^ 



x+y 
V2v 2 



V3 = V4 



4=f and fi = f 2 . 

s/2v f 



Identical calculations as given above lead to the expres- 
sions for Gfi(r,T) and Gt 2 (r, r). Putting all of them 
together we finally obtain 



G e (r,r) = 9(-r)e 



(H+Ah)\r\ m b 



G 6l (r,r) = ie(-r)e^ +4 ^)l-l x 



2tt|t| 

4irvfV 2 (t 2 + f f ) i 4,TrvfV2 (r 2 + f 2 .)! 



(38) 
if 2 ] 



/ 



d 2 fc 
(2^) 2 



Ak-r 



(b + ir)b 



-t X (ka) 2 \r\ 



e(-r)e^ +4 * x )l r l 



where m& 



2t X aS 



2tt|t| 

is the band mass of the boson. 



(34) 
(35) 
(36) 



iQ4-r 



4nVfV2 (r 2 + f 2 )2 4irvfV 2 (r 2 + f|) 2 I 

At this point let us make the following replacement in 
the spinon part of the physical hole propagator. 



1 



is nothing but the propagator for a particle of mass nib 
in 2d. Fourier transforming Gt^(k,U! n ) with respect to 
frequency yields 



(r 2 



(r 2 +r 2 Y 



(39) 



G n (k,r) 



1 



-E f \r\ 



sign(r) 



6f - lT)f 



E 



f 



This procedure is motivated by revisiting the following 
first quantized path integral formulation for the electron 
Green's function 



(37) 



G(r,r) = A" 1 / DaDx f Dx b e 



Next we expand the dispersion about the node 

Qi = (55^) 

e/(Qi+p) = 2xJa(p x + p y ) = Vfpi 
?7/(Qi+p) = -2AJa(-p x +p y ) = v 2 p2 

where Vf = 2\/2x Ja v 2 — 2^/2 A Ja and we resort to 
the notation p — p to denote 2d vectors. 



-iJ^- T} dr{L,+L b +L a ) 



i f, a-dx—i f, a-dx 

e 'b >f 

where integrates along the boson path X(,(t) and J 
integrates along the fermion path Xf(r). Lf^a are La- 
grangians for the fermion, the boson, and the gauge fields. 
Let I be the straight line connecting (0,0) and (r,r). 
Since the fermion is relativistic, the area spanned by the 



We now split up the momentum integral / &e** = fermion path \ } and | is of order L , where L 



2^2 



(27T) 

Iq (2J) 2 e t ^'" l "^' f where with i = 1...4 corresponds 
to (tt/2, tt/2), (-tt/2, -tt/2), (tt/2, -tt/2), (-tt/2, tt/2) in 
this order and we have set a the lattice spacing equal to 
unity. 

The momentum Fourier transform about Qi is then 
given by 

d2 P r t(Q*+p)-? 1 r -J(vfPi) 2 + (v2P2) 2 \r\ 

iQi-f r r . 



r +vir 

and f/ is the fermion velocity. For the non- relativistic 
boson, the area spanned by the boson path ^ and | 
is of order v^t^Jt Jmb, with mi, the boson mass. For 
t > \/m v 2 j, the area between the loop spanned by If 
and lb mainly comes from the area between ; f and | and 
we can approximate by Jj : 

G = N- 1 



sign(r) 



with cci 



^/(v/pi) 2 + (V2P2) 



- ix 1 - yi\ 



G fG 



4irVfV2 (r 2 + f 2 ) 2 



(40) 



Analogously we can obtain the contributions around 
the other three nodes - hence 



y— g 



G /T (r,r) 



rf2fc r ik-P l -E,\r\ 

(2tt) 2 2 



sign(r) 



ir) f 



Ef 



47rw/i;2 (r 2 +f 2 )2 
«2/3 + £3] 



e »Q 3 -f r r 



ATTVfV 2 (r 2 +f 2 )2- 
e iQ 4 -r j r + — 4 _ £ 4 j 



• f(r,T) 11 j ,j \ ^ f a-dx—i f a-dx 

A~ J / DaDx f e~^° dr(L f +L a ) e j, j (/ 



where Gf,o is the free boson Green's function (not cou- 
pled to the gauge fields). Note that the fermion Green's 
function G/ is just the gauge invariant Green's func- 
tion discussed in the first part of the paper. Hence we 
see that by performing the above replacement (|3^) with 
a = 27^ ~ 0.54 we make the transition from the free 
spinon Green's function to the particular gauge invariant 

Furthermore 



spinon amplitude discussed in section III 



4TTVfV 2 (r 2 +f|)2 4irv f v 2 (r 2 +r|)T 



we note here that our previous attempt of identifying the 
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FIG. 4: The top figure shows the good agreement between 
a direct evaluation of the convolution integral in momentum 
space (dashed line) and the result based on equation (^J) the 
Fourier Transform with a = 0. The bottom figure depicts 
the evolution of the spectrum as we increase a. The value 
a — 0.54 is "motivated" by our earlier analysis. Note that 
here t/J = 3 and we have used the implementation of the 
hypergeometric function given in "Numerical Recipes in C", 
W.H. Press et.al Eo 



physical hole Green's function with the gauge invariant 
spinon amplitude corresponds to neglecting the pj term 
in (|38|) which arises from the fluctuations of the holon 
about the classical path (the static straight line in the 
m& — > oo limit). In what follows we will analyze the full 
expression (^8j) under the replacement (^) where the 
fluctuations of the holons are assumed decoupled from 
the gauge field (hence the free boson pj contribution 
from the fluctuations around the classical path). This in 
particular implies that the area between the loop spanned 
by It and lb is solely given by the relativistic spinon path 
combined with the straight line segment which from the 
previous discussion is only exact in the limit m& — > oo - 
the static holon case. Thus the value a — 2^ ~ 0.54 re- 
flects the total gauge fluctuation effects strictly speaking 
only for the static holon case. Granted these limitations 
of the present analysis let us press on and analyze the 
effect of a finite positive a. Note also that we introduced 
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FIG. 5: Top: Direct evaluation of the convolution integral 
in momentum space (dashed line) compared with the result 
based on analytic continuation of equation with a — 0. 
Note that the result of equation ( jig ) actually goes negative for 
small uj which we believe is an artifact of the numerical eval- 
uation where we sum terms that alternate in sign but don't 
form a null sequence. This is supported by the agreement 
with the direct evaluation of the convolution integral which 
is small as can be seen in the figure in the region for uj — > 
where the numerics becomes tricky. The curve is obtained 
with n = 4 in the evaluation of (p8|) which shows the fast 
convergence for intermediate values of uj. The bottom figure 
depicts the evolution of the spectrum as we increase a. Note 

2 

that the onset of the spectrum is determined by 



r with units of energy to balance dimensions whose value 
will be on the order of the lattice scale. 
Thus our task now is to calculate 



G e (k,uj n ) = / cfrdre 



G e (f,r) (41) 



To do this we first concentrate on k = Qi +p and note 
that the other three points are related by symmetry. 
From expression ( |38| ) we see that only the first term in 
curly brackets contributes since all the other terms have 
large phases even in the infrared limit. We are thus left 
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with the evaluation of 
G e (Q x +p> n ) = J d 2 r dTe-^V^+l"'- 4 ^ 7 ' 



x e 

2ttt 



'V2« 



v 



4tt(t 2 v 2 + r 2 )- 



(42) 



where we have simplified to the case Vf = v-i — v and 
have taken care of the fact that our analysis only allows 
us to calculate the hole spectral function which was in- 
dicated by the 9(— r) factor in the boson propagator. 

Next let us analyze expression ( |4^ ) in some limiting 
cases where the integrals can be performed exactly. The 
general result and some more technical points can be 
found in the appendix. 

The simplest point to analyze is the expression right 
at the node p = where the terms proportional to xSzy 
average to zero on angular integration and we are left 
with 



oo />oo 



G e (Qi,o; n ) = K x I drr / (fre-l^+H- 4 ^ 
Jo Jo 



e 2t 



(t 2 v 2 + r 2 ) 



where K l = 
evaluated via 



L = v/T. The r integral can be 



Gradshteyn Ryzhik 5 edition (GR5) fig 3.382 (4) 



dx(x + pfe-^dx = -^e^T{y + 1, j3p) 



\arg/3\ < it, Rep, > 



where T(x, y) is an incomplete Gamma function. The 
resulting r integral can again be performed with the help 
of (GR5) 6.455 (1) 



r(f, ax)dx 



a v T(p + v) _ / 1 j8 

Re(a + [3) > 0, i?e^i > 0, Re{ji + u) > 



where 2 -Fi is a hypergeometric function and we have the 

~ 2 

following identifications (3 — [iw n + \p\ — 4tx\ — a = 

m b v 2 _ a+X _ a—X 

2 h 1 — 2 ~ 2 ' 

arrive at 



Putting this together we 



G e (Qi,w n ) 



m b r Q r(a) 



4tt (a + 1) [iw„ + \p\ - 4tx] a 



(43) 



a + 3 
2 fi l, a; — — ;1 



rribV 



2[iuj n + \fJt\- itx] 



We are now ready to perform the analytic continuation 
iu n — > lu — i0 + to obtain the hole spectral function at the 
node Qi 



a_(Qi,w) = — — - — — e(-uj) 

iir [a + 1) 
Re 2 Fx ( l,a;2±^ ; i 



sin(ira) 



[H-H + 4*x]" 0, e(H-H + 4tx) 



rribV 



2[u + \fi\ - 4ix - t'0+] 



(44) 



cos(7ra) 



M- | M | + 4< x r a ©(M- H + 4i X ) 



Q! + 3 

Ira 2 F\ \ 1, a; — - — ; 1 



2[u + |m| - 4t X - «)+] 
I 



Note that in principal we also have a contribution from 
0(|//| — 4tx +w) however since the chemical potential is 
exponentially close to the bottom of the band we will set 
\fi\ ~ 4tx for our calculation. It is nevertheless important 
in the derivation of the above that \p\ — 4i% = + for the 
integrals to converge. In Fig.Q we show the effect of a 
non-zero a on the spectrum at the node. The top figure 
simply compares the two results for the spectrum based 
on a direct evaluation of the convolution between single 



boson and spinon Green's functions (dashed line) and 
the evaluation of Eq. (E3) for the case a ~ where they 
should agree. The bottom portion of Fig.Q depicts the 
evolution of the spectrum as we increase a away from 
zero. Though still rather incoherent a shift of spectral 
weight toward lower energies is apparent. This piling 
up of low energy spectrum is directly associated with a 
nonzero alpha which in turn is a consequence of the gauge 
interaction between the holon and the spinon responsible 
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for the incoherent mean-field spectrum. Thus we see that 
when our reference system is the incoherent spectrum 
associated with the independent propagation of spinon 
and holon the gauge fluctuations have the desirable and 
intuitive effect of binding the two degrees of freedom back 
into the physical hole. This should be compared with our 
initial attempt (in the first half of the paper) to reference 
the gauge fluctuation effects with respect to simply a free 
spinon plus a coherent boson. A nearly coherent boson 
(with a well defined phase) is however not consistent with 
a massless gauge field whose origin lies in the strong local 
phase fluctuations of both holon and spinon. 

Fig. J5| shows a similar comparison for momentum p x = 
p y = —0.1 along the diagonal. This result is based on the 
general equation derived in Appendix C. The evolution 
with a follows the expected trend. The onset of the spec- 

2 

trum is governed by which gives rise to the peak in 
the spectrum centered at the boson energy. 



V. EVOLUTION WITH m b 

A finite a leads to a piling up of spectrum at the lower 
edge. However as we pointed out above this onset is gov- 
erned by the boson energy since the spinon spectrum is 
linear with velocity v the minimum of the spectrum for 
small momenta is always going to be determined by the 
boson with its quadratic dispersion. This of course re- 



flects our cavalier treatment of boson-boson correlations 
which in the absence of a full theory for the spinon holon 
+ gauge system we have totally neglected. Boson cor- 
relations will tend to reduce the density of low energy 
boson states which are responsible for the large distance 
behavior of the boson correlater. On the crudest level 
we might hope to capture this reduction by simply de- 
creasing rat which determines the free boson density of 
states. 

In order to see this dependence most clearly let us re- 
derive an expression for the spectrum at the node. In- 
troducing t — tv we would like to evaluate 



G e (Qx,cj n ) 



df / d 2 re-^+^- ii ^ 
i] Jo 



8tt 2 



\f 2 +r 2 ] 



We can now view this as a 3d integral over the positive 
hemisphere where f = RcosO r = RsinO and after inte- 
grating out <j) and R we are left with 



— - / dx- 

v Jo [f(iw„ 



G°(Qi,w») 
T(a) 



After analytic continuation we arrive at the following ex- 
pression for the hole spectral function 



J 



A-( Ql ,c,) = e(-a,)^ Sfmrar(a) /'V- 



mbV 



2x 2 



[x 2 (|a,| - | M | + 4t x ) - - ^ )]c 



which in the a — > + limit reduces to 



VI. m b = 



e(-w) 



-dx 



lim 

w-*||t|-4tx 



2[M-| /i |+4t x + m b „:* / 2] 

K x |w|-|M|+4t X 



(45) 



We thus find that — determines the slope of the spec- 
trum at the onset. In Fig.|7[ we depict the effect of de- 
creasing the mass of the exponential by a factor of 3 
(t = 9 J) and 4 respectively for the case a = 0. Note that 
in addition to becoming steeper in accord with the result 
just derived for the case p = the onset of the spectrum 
is pushed to higher frequencies initially simply following 

2 

the increase in The curves seem to converge to a 

step function, however, the onset of the spectrum does 
not converge to pv the energy of the spinon but rather 

pv 
~2~- 



As was indicated in the previous section the free boson 
dispersion determines both the onset an the shape of the 
low energy spectrum. Correlations will however reduce 
the low energy modes whose effect we tried to emulate 
by reducing the boson density of states. Let us now set 
To;, = from the start. The rational behind this section 
is the attempt to capture the correlations corresponding 
to collective bosons that are not condensed. For nearly 
condensed bosons the appropriate collective variables are 
the density and phase and we shall assume that we can 
neglect fluctuations in the density and only the correla- 
tions of the phase field determine the infrared physics. 



(&t(r,r>(0,0)) 



where the dynamics of 9 is given by a 2d AY-model. 
The 8 Green's function can thus be determined and the 
resulting boson correlation function reads 



(6 f (r, r)b(0, 0)) = p e 4 ™->« ! e 4 ' VbXi V r 



(46) 
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FIG. 6: Linear approximation Eq.(^5|) compared with full 
result for the case a = 0. 



which in the infrared limit we can write suggestively as 

1 \ 



(&t( T ,r)&(0,0)) ~p e *™ bXbl 1 + 



Xb denotes the boson compressibility Vb is the velocity 
of sound of the collective bosons and I is some lattice 
scale cut-off determining the wavefunction renormaliza- 
tion which will give rise to a reduction of the superfluid 
density pq. 
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FIG. 7: Evolution of spectrum with mi, at p = (—0.1,-0.1) 
at a = 



(47) 



We will take the second term in parenthesis as reflect- 
ing the uncondensed but correlated boson contribution 
to the long distance behavior of the boson propagator. 



Taking the above correlation function as the holon part 
of the physical hole Green's function we are faced with 
the evaluation of 



G e (Qi+p>„) =L~ 



d l r I dre-^e-^+^-^poe'^n^ 



; x±y 
'V2v 



^KVbXb\Jr 2 + v 2 t 2 Att(t 2 v 2 + r 2 )' 



r 



Let us first use the exponential parametrization 
1 1 f°° 



(r 2 +v 2 T 2 Y T{p) J 

which allows us to perform both the r integration via 
(GR5) 6.631 (4) 



drr v+1 J v {(3r)e 



ft 2 



and the r integral via (GR5) 3.462 (1) 



POO / 

/ x^e'^-^dx = {2l3)- y ' 2 T{v)exp[ 



x D. 



where D v is a parabolic cylinder function. 



Once the r and r integrals are done we are left with the 
integrals over the parameters s, u. Concentrating for the 
moment on the term proportional to r in the numerator 
we obtain after r and r integration 



2irK 



r(^)r(|)7 

2 [i^„ + |n|-4tx] 



dsdus x l 2 u^ 



e 4 <=+") e 



-Dn 



ImI - 4tx] 



4(s + u)(sv 2 + uv 2 ) " 2 



J b t" u " 7 > \ y/^( sv b + 111,2 

Next we perform a change of variables 

v 2 s = t\Si v 2 u = (1 — ti)si 
which changes the parameter integration to 



2-kK 



r(^)r(|) 



dh 



2„2 



V V 



-1/2 
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xexp 



'(l-ti)si 



[iuj n + - ttxf 
8si 



4(tl/«f + (l-tl)/« 2 ) = l 



4s 



2 / ti , 1-ti 



/2sl 



Finally we can perform the integral over s\ via (GR) 
7.725 (6) 



D- V [2(kt)*]dt 



2^-fV^(/3) i (v /3. „ + /?+l .z-fc 



2 1 2 

with the identification [3 



2 fc 



, + |M|-4t X ] 2 
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turns 



out that for the numerical evaluation of the final param- 
eter integral over t\ it is easier if we use the following 
quadratic transformation formula for the hypergeomet- 
ric function 



F(a, b; a - b + 1; z) = (1 - zy a x 
fa a , 1 , 4z 
F U'2- 6+ 2 ;a - fe+1; -(I-^ 



Putting everything together we obtain for the hole 
Green's function 



G e (Q l +p,oj n ) = K{ 



dtt~^(i - 1)~ 



tv 2 + (l - t)v£ 



2 \ ~ 



V 2 V 2 



[iu) n + \fi\ - 4tx] 



tv 2 + (l-t)v 



2 F\ 



( a _ 1 3 + a [iu, n + \»\-4t X ] 

ry • • 



-a-«K p \ 

av(p x + p v )V2 



V 



tu 2 +(i-t)ti 2 



dtr*(i - 1)- 



tv 2 + (1 - i)^ 



2\ ~ 



fu 2 + (l-t)tj 2 



1 + Q 



2 ^1 



1 + a 



a + 1 3 + a 



[iw n + |mI - 4 *X] 



to 2 + (l-t)-u 2 



!!-■('- 



A) 



2 ' 



[iw„ + |mI - 4*x] 

r 



to 2 +(i-t)-» 2 



(48) 



and if 



A°p e 47r "t.x b ' r(q) 



. We can now continue 



u — i0 + analytically and take the imaginary part 
to obtain the spectral function. In Fig.[|| and Fig.|| we 
plot the spectra for two momenta at p x = —p y = —0.1 
and p x = Py = —0.1 at different values of a. As opposed 
to the spectra with finite boson mass Fig.|| and Fig.||] 

show the spectral onset at the spinon energy v*/p 2 + p 2 



We can easily understand this by noting that v 2 = — 
which with the mean-field values for /i and is equal to 

2 



This in turn implies that for a given momentum 

the onset in the spectrum occurs when all the momentum 
is taken up by the spinon resulting in a lower edge of 
the spectrum tracing out the spinon dispersion |^o| . As 



with the case of finite boson mass the spectra for zero 
a are broad without any special features at the onset 
of the spectrum. The confining tendency of the gauge 
fluctuations can again be nicely observed on increasing 
a. 



VII. CONCLUSION 

The topic of this paper - the single particle spectral 
function - was analyzed in a myriad of ways. After briefly 
addressing the problem of gauge invariance in the con- 
text of the single particle Green's function we gave a first 
quantized analysis of a particle coupled to gauge fluctu- 
ations. The proposed gauge invariant quantity to reflect 
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FIG. 8: Evolution of spectrum with a at p x — —p y = —0.1. 
In order to get reasonably smooth graphs we first calculated 
the spectra with a resolution Auj = 0.001J and then averaged 
5 points to plot the spectrum respectively. 
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FIG. 9: Evolution of spectrum with a at p x — p y — —0.1 
along the diagonal. We averaged 10 data points to plot the 
spectrum which corresponds to Aui = 0.01J 



the propagation of the particle was a loop, one segment 
of which was made up by the particle path propagating 
from fi to ff ■ The loop was completed by a straight line 
segment connecting the two destinations. We pointed out 
that if a particular gauge (an axial gauge with a ■ n = 
with ft oc ff — fi) was chosen for the gauge propagator 
the contribution of the straight line segment to the par- 
ticle amplitude could be made to vanish. With this idea 
we then analyzed the amplitude in second quantized for- 
mulation and obtained the dressed Green's function with 
a slower decay in space-time. However this approach has 
some spurious infrared problems related to the incom- 
plete gauge fixing in the temporal gauge and required 
special (sometimes ad hoc) regularizations (Lcibbrandt 



prescription). To confirm the result, we then resorted to 
discuss the above amplitude as a gauge invariant two- 
body propagator with one "body" infinitely heavy and 
hence giving rise to the straight line path introduced by 
hand earlier. Thus freed from choosing a particular gauge 
we performed a proper analysis of this propagator. The 
resulting anomalous dimension for the gauge invariant 
amplitude agrees with the previous calculation. It how- 
ever could not and should not be interpreted as repre- 
senting the propagation of a physical particle. 

This forced on us a reinterpretation of the gauge field 
- not primarily as a scattering mechanism - but rather 
a confining field which wants to form a bound state be- 
tween the two charged particles propagating under the 
mutual attraction mediated by it. This change in per- 
spective led to the real space analysis of the second 
part where we incorporated the gauge fluctuation effect 
heuristically into a slower decay of the spinon contribu- 
tion to the hole propagator. We saw how this change in 
the spinon Green's function led to the piling up of low 
energy spectral weight. We associate this with the con- 
fining tendencies of the gauge field. Our treatment of the 
bosonic holons as free, however clearly overestimated the 
density of low energy holon states which led to the spec- 



trum to be peaked about the holon energy 



This is 



not seen in experiments ( with the possible exception of 
the single hole spectrum which seems to have a parabolic 
band near (n/2, 7t/2) jSlJ compared with the Dirac spec- 
trum observed at higher doping levels) . In order to over- 
come this we argued for the holon Green's function to be 
determined in the infrared limit by collective correlations 
rather than the free particle behavior. This remedied the 

2 

problem of the edge of the spectrum tracing out and 
led to the spectra depicted in Fig.Q and Fig.||. It is 
however clear that at the level of our present study the 
real value of a is not known and hence we are forced to 
treat it as a phenomenological parameter. A proper mi- 
croscopic analysis of the holon + spinon + gauge system 
is needed to put the heuristic results proposed in this 
paper on a firm foundation. 

We would like to remark that the electron spectral 
function discussed here is the high energy spectral func- 
tion, which is valid for high energies and/or high tem- 
peratures. This is because we have ignored the boson 
condensation and the instanton effect. For T < 150K, 
the boson-condensation/instanton-effect will be impor- 
tant which will lead to a spin-charge recombination |)| 
due to the confinement of U(l) gauge field. In this case 
a 5-function peak will appear in the electron spectral 
function, and low energy excitations will be described 
by electron-like quasiparticles. However, the Fermi sur- 
faces of those well defined quasiparticles are formed by 
four small segments. [|] 

We would like to thank P.A. Lee, T. Senthil and A. 
Seidel for helpful discussions. This research is supported 
by NSF Grant No. DMR-01-23156 and by NSF-MRSEC 
Grant No. DMR-98-08941 
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VIII. APPENDIX A: FERMION SELF ENERGY IN THE a = GAUGE 

In this appendix we present details of the evaluation of the self energy in the ao — gauge at T=0 



N 



d 3 q 



(g -Pa)q 2 + 2(q -p )q 2 



cpip- q) 2 ql 



d 3 q f (p x - q x )(ql - ql) + 2q x q y (p y - q y ) 



(2tt) 3 



d 3 q f (Py - <ly)(Qy ~ ll) + ^QxQyiPx - Qx) 



(2tt) 3 



(49) 



Let us take a closer look at the 7°-term which has two contributions 



d 3 q 
(2tt) 3 



(go - Po)f 

2 n 2 



'(p(p- q) 2 q 2 



(go - Po) 

/¥(p- q) 2 



(50) 



First look at term B which can be calculated in the standard way using dimensional regularization combined with the 
Feynman trick in the generalized form 



1 



z a - l (l- z) - 1 T(a + P) 
A a BP ~ J Q [zA + (1 - z)B] a +P T{a)T{(3) 



dz 



With A = (p- q) 2 B = q 2 we find for 



d 3 q 2(q - p ) 
(2tt) 3 ^/f(p- 



| div 



Po 



6tt 2 3 - d 



Term A exhibits in addition to the UV divergence which presumably will be taken care of by some appropriate 
regularization procedure an infrared divergence due to the q 2 , in the denominator. This divergence is an indication that 
our gauge fixing ao = has a large residual gauge freedom in the form of time independent gauge transformations. 
This reflects the fact that in using the temporal gauge we have "lost" Gauss' law which within the Lagrangian 
formulation is obtained as a constrained on the gauge fields by varying the action with respect to ao. In order to deal 
with this problem we follow Leibbrandt and Staley p2]| which have introduced a prescription for this unphysical pole 
in the form 



1 



where the limit fi — > is taken after the integrations over loop momenta and Feynman parameters is completed. This 
leaves us with the following expression for A 



= I 



d 3 q (go - Po)f 



(2irY ^f(p-q) 2 {qt+V 2 ) 2 
To compute this integral (23) we first combine the denominator using exponential parametrization 

i r°° * 3 



(51) 



(p- q) 2 V<P(Qo + M 2 ) 2 



dtidt2dty, 



-.exp [-h(p - q) 2 - t 2 (f - t 3 (q 2 + n 2 )] 



then we use dimensional regularization to compute 



/ 



d d q 
{2n) d 

1 



d 



d(h+t 2 ) dt 3 

d/2 



d(h+t 2 )dt s (27r) d \ti + 1 2 J {h+h+h) 1 / 



3 [-(ti+t2)9 a +2tii?g-t 3 <Zo] 



(*l+*2) (*l+*2)(*l+*2 + *3) 
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In order to compute the divergent part of this integral we concentrate on the terms with no powers of p in the 
numerator from the differentiation since it is not hard to convince oneself that all terms of higher order in p converge 
in d = 3. Thus to extract the divergence we have 

J (27T) 



Po 



q % 



O 3 V¥(p-q) 2 (qi + v 2 ) 2 

v^oU « 2 ^ 



1-d 
2 



(h+t2)- (1 - d)/2 3 {t 1+ t 2 )^- d )l 2 



(tl + *2 + i 3 ) 3 / 2 2 (t! + t 2 + t 3 ) 5 / 2 



tip 2 



ht 2 P 2 



hp - p h 



Changing variables to 
we have 



(tl+t 2 ) (tl+t 2 )(tl+t2 + ts) 

ti = siu t 2 = st(l — u) < 3 = s(l — t) 

^L^l j 1 du f'df [ X ds a 3-l/2 e - a [ t (l-«)„^(l-t)t« a P S + ,. a (l-t)] 

2^ (2n) d J J J 



1-d 



-(2+d/2) 



(1 - i)i" d / 2 3 



The integral over s yields 



T d/2 



l>{) ::' " (3-d 



2V^ (27r) d 
3 



(1 - u)V2 2 

1-d 



- _ s -(2+d/2) f l-d/2. 



(1-i) 



1 ,1 

dw / dt 
(I Jo 



(1-U)V2 

(l-tjt-^a-u)- 1 / 2 - 



t l-d/ 2 (l_ t )(l_ u )-l 



/2 



{t(l - ti)wp 2 + (1 - t)tu 2 pl + /i 2 (l - t)} 



(d-3)/2 



Expanding the term {• • -}( d 3 )/ 2 for d — > 3 and integrating over the remaining two parameters we finally obtain 



y (2tt 



In a similar fashion we arrive at 



and thus finally 



d 3 q q a c? 



\div 



Po i_ 

2tt 2 3 - d 

Pa 1 



(^) 3 V¥(p-q) 2 q 2 ^ 3 ~ d 



N 



i^o\div = 70 



3p 1 Po 1 



2tt 2 3 - d 6tt 2 3 - d 
■• .. ' .. ' 



= 3^ 7 P °—d 

The spatial components of the self energy can be evaluated similarly to give 



N 



2— x,y 



(52) 
(53) 



where the following results were used 



d3q 


, 




1 


(2tt)3 , 




2tt 2 3 


-d 


d 3 q 




Pi 


1 


(2tt)3 , 




3tt 2 3 


-d 


d 3 Q 


^ 1 - 


1 


1 


(2tt)3 , 


/¥(p-q) 2 q 2 


2vr 2 3 


-d 



with j, fc = 1, 2 the spatial components 
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IX. APPENDIX B: STATIC BOSON SELF ENERGY 

To calculate the static boson self energy we need to evaluate 

-T, b {p )= / 7^2—7 r- (54) 

J (2tt) 2 t(p Q - g ) 

where Doo = j^-^r in the Landau gauge. To combine the denominators we use the following variant of Feynman's 
trick 

00 „f ) - 1 r 



a a bP r(a)r(/3) J {a + yb) a +P 

which in our case gives 

dy 



\q] 3 ( Po -qo) 2j (<?+y(p -qo))V 2 

Combining the denominator (q — y/2v) 2 — y 2 /4 + ypo with v — (1, 0, 0) and shifting the momentum to k = q — y/2v 
we obtain 

N . . 3 f°° f d d k k 2 

1— T, b (po) = - dy jT—rj-r; — 7 

which with 

d^k kpk v 1 8^ v L( 2 ) A L 



{2*) d [P + A]V2 2(4tt)3/ 2 r(|) 



A" 



y 2 

A = ypo - - j 
The final parameter integral is calculated using (lj] 

y (ay + 6 )^=_ _ 

and yields 

y 2 sd-i PoT(3-d) 



^ dy(ypo - —) 



3 



r(^) 



4 v 2 



which finally gives 



X. APPENDIX C: FT OF THE HOLE GREEN'S FUNCTION 

In this appendix we analyze 



G e (Qi+p>„) = I d 2 r I d T e-^-V^"+l^- 4 ^ T ^— e" 1 ^^ V ' 2v 1 3 _ a (56) 



V2v 



lo 2-kt Ak(t 2 v 2 + r 2 f 

more fully than was done in the main body of this paper. Let us split the above integral into three parts 

P POO T — a 1 

h = d 2 r dTe -W-r\-[^ + M^i X ]r^_ e -^- 



o 2-kt 4tt(t 2 v 2 + r 2 ) 

/pOO t — ct 2 — % — 
d 2 r / dre-^e-^+^-^^^e-^ 
Jo 2-kt 4 7 r(r 2 i; 2 +r 2 ) i 

Is = h(x^y) 
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e 2t 



and first concentrate on Ii After angular integration we arrive at 

/>oo />oo 

h = Kj dr drre-^+M-^ J (pr) — 

Jo Jo (t z v z + r z ) 2 

with J the zeroth order Bessel function and K\ = rrLbV , I l — . Next we utilize the series representation of Jq 



J v (x) = ( -x\ 



k=0 



k\r{v + k + i) 



which converges absolutely for all x. This fact allows us to integrate the series term by term 

1*00 

dTe -[iv n + M-4txlT( TV ja-3 / dxx k e -^ 



00 (-ip 2 ) fc 

Using GR5 (3.383 (5)) 

POO 

/ e' vx x q - 1 {\ + ax)-"dx = a- q T{q)^{q,q + 1 - v;p/a) Req > Rep > Rea > 
Jo 

where \I>(a, b; z) denotes a confluent Hypergeometric series. The r integration can again be performed with the help 
of GR5 (7.621 (6)) 



/•OO 

Jo 



e _ stdt = r(b)T(b c + l) _ b 2Fl ( a b . a + b _ c + 1 . 1 _ s - 1) 

I (a + b — c + 1) 
Reb > Rec < Reb+1 



and the identifications b = 2k + a a = k + 1 c = k + i±« s = 2[iu n + — At~x\/ {m^v 2 ) yielding to 



fe=0 



r(2fc + a)r(fc + i±g) 
r(2fc + 3P) 



.,r,(fc + l,2fc + a;2fc+^;l- ™*> 

2 2[iLu n + \n\ - 4tx\ 



(57) 



Notice that this result reduces to the one quoted in the main body of the thesis for the case p = where only the 
k = term contributes in the sum. 

Next let's look at I 2 which can be cast into 



I 2 



TOftL 



dr d 



AV2n Jo t dp x J 



With -^-J (p\/x) = — J\{p\fx)p x jp the above program can now be repeated step by step. Thus combining Ii,l2 and 
I3 we find for the physical hole Green's function 



G e (Q!+£0 = 

- mbTa £ + imi - 4*]— T(2k + a)T(k + m 



(58) 



8tt ^ k\ 

k=0 



r(2fc + 2±«) 



2 F (fc + l,2fc + a;2fc+^;l- "*» 



8 - V2 ^ A ' ! 



E 



/ 2 2 \ 



5 + a 



_ i r(2fc + a + l)r(fc+i±S) 



r(2fc + ^s) 



2F1 fc + 2,2fc + a + l;2fc + 



mi,v 



2[iu n + |m| - 4t X ] 
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